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ABSTRACT 

We present narrow-band images of the Red Rectangle (RR) nebula which reveal the distinct morphologies of 
this intriguing nebula in different optical emission bands. The morphology of the RR nebula in blue lumi- 
nescence (BL) and extended red emission (ERE) are almost mutually exclusive. We also present the optical 
detection of the circum-binary disk of the RR in the light of the BL. The total intensities from the two optical 
band emissions (BL and ERE) when summed over the nebula are of comparable magnitude. Their spatial dis- 
tributions with respect to the embedded illumination sources lead us to suggest that they may be attributed to 
different ionization stages of the same family of carriers. 

Subject headings: ISM: individual(Red Rectangle) — ISM: molecules — dust, extinction — radiation mecha- 
nisms: general — stars: individual(HD 44179) 



1. INTRODUCTION 

The optical spectrum of the Red Rectangle (RR) neb- 
ula consists of several components of comparable intensities 
(Schmidt et al. 1980; Witt & Boroson 1990; Vijh et al. 2005), 
and it is possible to use narrow-band imaging to isolate the 
morphological substructures where particular emissions are 
dominant. This is in contrast to most ordinary reflection 
nebulae, in which the scattered light component dominates 
(Witt & Boroson 1990). The special status of the RR is a con- 
sequence of the optically-thick circumstellar disk surrounding 
the central illuminating source and seen edge-on, which pre- 
vents forward-directed scattered radiation from reaching the 
observer. Thus, it is the suppression of scattered radiation by 
at least one order of magnitude that makes blue luminescence 
(BL) and extended red emission (ERE) readily observable in 
the RR. The latter two radiation processes are most likely 
the result of isotropic emission processes, which are readily 
observable from nebular locations where the scattered light 
is preferentially directed away from the observer on account 
of the strongly forward-throwing scattering phase function of 
dust at optical wavelengths. 

The spectral components of the optical light of the RR are: 
a) Dust- scattered radiation: to first order, this spectrum is a 
somewhat bluer version of that of the more luminous star in 
the central binary system, HD 44179, a T eff -8250 K A-giant 
AGB star, b) Blue luminescence (BL) band: this emission 
band with a peak near 380 nm (FWHM ~ 45 nm) has been 
attributed to fluorescence by neutral PAH molecules with 14 
- 18 C atoms. (Vijh et al. 2004, 2005) c) Extended red emis- 
sion (ERE): this band peaks near 670 nm (FWHM - 180 nm) 
and it is the primary cause for the unusual red color of the 
RR. Although observational constraints have indicated that 
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the ERE carrier requires photons with E > 10.5 eV for their 
initiation, indicating a possible association with doubly ion- 
ized molecular species (Witt et al. 2006), the carrier of the 
ERE is as yet unidentified (Witt & Vijh 2004). d) Sharp red 
emission features: this is a set of relatively narrow emission 
bands due to as yet unidentified molecules. They appear in 
a wavelength range largely coincident with that of the ERE 
band (Van Winckel et al. 2002). It has also been suggested 
that some of these bands are emission counterparts of cor- 
responding diffuse interstellar bands (DIBs) (Scarrott et al. 
1992). Spectroscopic studies have demonstrated that the in- 
tensities of the emission bands b), c), and d) are comparable 
to those of the underlying scattered light continuum, allowing 
their relative distributions to be traced. This makes it pos- 
sible to constrain the physical state(s) of the respective emit- 
ters, because the almost-edge-on geometry of the bi-polar out- 
flow allows one to associate specific physical conditions with 
regard to density, radiation field, and ionization equilibrium 
with specific nebular regions. 

Our imaging study of the RR differs in significant as- 
pects from previous efforts. Most earlier work was fo- 
cused on sub-arcsecond resolution at near- and mid-infrared 
wavelengths with the aim to resolve the core of the RR 
(Roddier et al. 1995; Cruzalbes et al. 1996; Mekarnia et al. 
1998; Tuthill et al. 2002; Miyata et al. 2004), which resulted 
in the well-resolved imaging of the circumstellar disk and 
brightness distribution of the RR in the bright innermost re- 
gion. The more extended brightness distribution of the RR 
was studied at high resolution by Cohen et al. (2004), using 
WFPC2 on HST. Their band-passes were limited by the fil- 
ters available in WFPC2, with the F467M filter defining their 
shortest-wavelength band. Thus, they were unable to ad- 
dress the spatial distribution of the BL, which peaks below 
400 nm in wavelength. They were, however, able to demon- 
strate the distinctly different spatial distribution of the dust- 
scattered radiation in the RR, which dominates the F467M 
band, and of the ERE, which was imaged in three different 
red bands (F588N, F622W, F631N). They also produced the 
best-resolved record of the ladder structure in the ERE cross- 
ing the outflow cones of the RR. Another important mapping 
effort is that of Bregman et al. (1993), who imaged the ex- 
tended brightness distribution of the RR in the two aromatic 
emission features at 3.3 fim and 11.3 jim. They demonstrated 
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TABLE 1 

Narrow-band observations of the RR nebula. 



Filter Filter width (A) 


Exposure (s) Dominant Component 


3934 27 


5x900 BLband 


4050 57 


5x120, 1x180 Scattered 


5700 135 


5x 120 Scattered + ERE 


6400 135 


4 x 60, 1 x 1 80 ERE + Scattered 



that the 1 1 .3 fim brightness distribution followed the extended 
bipolar structure of the RR, while the 3.3 /im isophotes were 
more concentrated near the center of the nebula. Our goal was 
to pursue the large-scale brightness distribution in the RR of 
specific optical band emissions, which we attempted to iso- 
late with narrow-band filters and still narrower Fabry-Perot 
bandpasses, while the angular resolution was limited by see- 
ing (~ 1 Our work includes, therefore, a first effort to 
image the BL in the RR. 

2. OBSERVATIONS AND REDUCTIONS 

The data were obtained at the Apache Point Observatory on 
January 7- 11, 2005, using the Goddard Fabry Perot (FP) on 
the 3.5-m telescope. Images of the RR were obtained using 
narrow-band filters centered at 3934 A ( FWHM ~ 27 A), 
4050 A (FWHM - 57 A), 5700 A (FWHM - 135 A), 6400 A 
(FWHM ~ 135 A). The nights were partially cloudy and only 
relative fluxes could be obtained from these images. We used 
the FP etalon to image the narrow emission features (item d 
of the list in the introduction), but could not get good enough 
exposures due to clouds. This paper deals only with the pre- 
pass filter images taken during the run. 

Data reduction was carried out using standard IRAF tasks. 
All images were trimmed to a size of 561 x 561 pixels or 
3. '4 x 3. / 4. The darks, twilight flats and object frames were 
corrected using an average bias frame. The flats and ob- 
ject frames were then dark-corrected using a combined dark 
frame. The flat frames for each filter were obtained by moving 
30" between successive exposures to avoid coherent addition 
of bright stars. These flat frames were then median-combined 
for each filter. Dark corrected object frames were then flat- 
fielded using the combined twilight flat for each filter. Object 
exposures were taken in sets of five, in a pentagram pattern of 
pointings. Object frames with the best seeing were selected 
and median-combined using tasks imalign and imcombine. 
Table 1 shows the exposure times for each filter. Stars in the 
frames were used to register and shift the images: this worked 
for all filters except the 3934 A filter. As no stars other than 
the central extended source were visible on these frames, the 
center of the extended source was used to shift and co-add 
the images. This was found to be reliable after the fact as 
one of the next-brightest stars in the field now appeared in the 
co-added image, albeit with extremely low counts. All im- 
ages were further cropped to 200 x 200 pixel or 73" x 73" to 
best show the RR. In the absence of calibration, images were 
scaled using background counts at large offsets from the cen- 
tral bright source before computing the ratio images discussed 
in the next section. 

3. RESULTS 



3.1. Blue Luminescence in the Red Rectangle 

The two narrow-band images at 3934 A and 4050 A (Fig- 
ure 1 and 2) show the appearance of the RR in the wave- 
length region of the BL band. The former is much closer 
to the peak of the BL emission, while the latter covers the 
long-wavelength tail of the BL band. This causes the A4050 
image to include a larger contribution, in relative terms, from 
the scattered light component of the spectrum. Figure 1 is an 
image of the RR taken with the 3934 filter, superimposed with 
intensity contours. The most remarkable feature of this image 
is the fact that the inner isophotes are consistently elongated 
in the E-W direction. A check on stellar images in the 3934 A 
frames confirmed that this elongation is real and not the result 
of guiding errors in the individual 900 s exposures. This E-W 
elongation in the blue image was confirmed in a subsequent 
observation obtained at APO on December 16, 2005. As these 
observations too, are not calibrated and less deep than those 
discussed in this paper, they have not been used in the anal- 
ysis. The E-W elongation is further illustrated in Figure 3, 
where we compare N-S and E-W traces of the intensity distri- 
butions at 3934 A and at 4050 A with the point-spread func- 
tion (psf) of the detector, derived from stellar images in the 
field imaged in the 4050 A band. In the N-S direction the in- 
ner core of the 3934 A image, which is more strongly affected 
by the BL than the 4050 A image, is barely resolved, in con- 
trast to the broader 4050 A image. This demonstrates that the 
BL is exhibiting a strong intensity peak at the position of the 
circumstellar disk, with a steep intensity gradient toward the 
regions above and below the disk. The 4050 A N-S profile, on 
the other hand, exhibits a comparatively broader distribution 
in this direction. A quite different result emerges from the 
traces in E-W direction. Both BL-dominated (3934 A) and 
scattered light-dominated (4050 A) intensity profiles of the 
RR are well resolved, and the ratio of BL to scattered light 
is increasing with offset from the center in both E and W di- 
rections. This set of traces confirms this same trend observed 
in the RR with long-slit spectroscopy with E-W oriented slits 
by Vijh et al. (2005), which showed that the BL emission ex- 
hibits a strong preference for the parts of the RR that are lo- 
cated in the disk, shielded from direct irradiation by the star. 

The E-W elongation of the BL disk is consistent with the 
identification of the BL as originating from fluorescence by 
small (3- to 4-ring) neutral PAH molecules as proposed by 
Vijh et al. (2005). Such molecules are photo-excited by stel- 
lar photons of ~ 4 eV energy, which can diffuse into the outer 
disk regions via scattering. The same molecules would be 
photo-ionized by 7.5 eV photons, which are present in regions 
of the RR with direct lines-of- sight to the central sources, e.g. 
the outflow cones. Once ionized, the same molecules assume 
a different electronic energy level structure and are no longer 
able to fluoresce in the blue spectral range. Vijh et al. (2005) 
examined other potential carriers of the BL band such as sil- 
icon nanoparticles, hydrogenated amorphous carbon and sil- 
icon carbide, and found their spectra to be inconsistent with 
the obsereved BL spectral characteristics. We believe, there- 
fore, that the 3934 A image is showing us emission from the 
outer portions of the optically thick disk surrounding the cen- 
tral stars in the RR, i.e regions where the local radiation field 
causes the ionization equilibrium to be preferentially neutral. 
We note that the isophotes of the 4050 A image, shown in 
Figure 2 are distinctly circular in the inner region of the im- 
age, in contrast to the elongated 3934 A isophotes. This re- 
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fleets the fact that dust- scattered light is dominant here. The 
outer isophotes of the 4050 A image become almost square, 
with the two sides being parallel to the polar axis of the RR 
nebula. While the BL image at 3934 A shows slight inten- 
sity enhancements at the base of the outflow walls on the 
W-side of the nebula, no such enhancements are apparent in 
the 4050 A image, despite the fact that this is a deeper expo- 
sure. The dimensions of the RR disk implied by the narrow- 
band imaging observations are in agreement with the results 
of Bujarrabal et al. (2003, 2005) who mapped this disk in CO 
emission with the IRAM interferometer array. Their map 
shows the disk having an angular extent of ~ 6". The 3934 A 
data show elongated isophotes over an extent of 6.5", which 
is in excellent agreement, taking into account the ~ \" seeing 
disk of the optical data. 

A perplexing aspect of the 3934 A image is the orientation 
of the disk (PA ~ 90°) compared to the orientation of the polar 
axis of the RR (PA ~ 11°). If the disk were exactly perpendic- 
ular to the polar axis, a PA of 101° would have been expected. 
This PA of 101° appears to be in agreement with disk images 
obtained in absorption at near-IR (Tuthill et al. 2002). How- 
ever, these near-IR images show the disk only over an extent 
of 200 mas to 500 mas, while the elongated emission recorded 
in the 3934 A image extends over a much larger angular ex- 
tent, larger by a factor of 10 to 20. The CO observations 
(Bujarrabal et al. 2005) also indicate that the CO disk is per- 
pendicular to the RR axis and is aligned with the disk orienta- 
tion derived from the near-IR and optical images. To establish 
that we would indeed be able to see a 1 1 ° rotation in our im- 
age we made a high-res test bar (0.3" high, 4" long, PA 101°) 
to simulate the RR disk. We degraded this test-image to the 
observed resolution by using the observed psf (at 4050 A) and 
the scattered light profile of the RR from the 4050 A image. 
The contours of the simulated bar clearly run at PA 101° and 
not in the pixel direction. It appears likely, therefore, that we 
are recording emission from the outer reaches of the circum- 



stellar disk, with the intensity distribution strongly influenced 
by the excitation conditions for the molecules responsible for 
the BL. 

3.2. Extended Red Emission 
The 6400 A image is dominated by the ERE (see Fig. 4). 
When compared to images at wavelengths < 5000 A, this im- 
age demonstrates that the X- shaped structures (variably re- 
ferred to as "arms" or "whiskers" in the RR literature) are 
a result of the ERE band. We interpret these structures as 
the walls of the internally illuminated bi-polar outflow cavi- 
ties, seen in projection. When divided by the 5700 A "con- 
tinuum" image, the resulting 6400/5700 ratio image, Figure 5 
(left) convincingly demonstrates that the ERE is preferentially 
emitted by material in and within the walls of the outflow cav- 
ities, and is extremely weak outside these regions. This is 
demonstrated even more clearly in the ratio image 6400/4050, 
Figure 5 (right). The 5700 A image is slightly overexposed, 
leading to bleeding from the stellar image and to the appear- 
ance of diffraction spikes that are not visible on the other im- 
ages. Figure 5 is basically a display of the ratio of ERE to 
dust- scattered radiation in the RR. We can interpret the rela- 
tive confinement of the ERE to the outflow cavities and their 
walls in one of two ways: either the carrier of the ERE ex- 
ists only in the outflow cavities and its walls, or the excitation 
requirements for the ERE restrict this emission process to re- 
gions which can receive the photons required for ERE exci- 
tation from the central sources. Either interpretation is con- 
sistent with the findings of Cohen et al. (2004), who demon- 
strated that the characteristic X- shaped structure of the RR is 
visible in several red bands dominated by the ERE but is ab- 
sent in the F467M image. 

4. DISCUSSION 

4. 1 . Excitation of the ERE 

The two possible interpretations of the apparent confine- 
ment of the ERE to the walls of the outflow cones of the RR 
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Fig. 3. — Intensity traces with the point-spread function (psf) of the detector. N-S (left) and E-W (right) traces of the intensity distributions at 3934 A and at 
4050 A are compared with the psf of the detector, derived from stellar images in the field imaged in the 4050 A band. In the N-S direction the inner core of the 
3934 A image, which is more strongly affected by the BL than the 4050 Aimage, is barely resolved, in contrast to the broader 4050 A image which hints at a 
comparatively broader distribution. In the E-W direction, both BL-dominated (3934 A) and scattered light-dominated (4050 A) intensity profiles of the RR are 
well resolved, and the ratio of BL to scattered light is increasing with offset from the center in both E and W directions. 




Fig. 4. — RR image at 6400 A, with 135 A wide filter. The image at this 
wavelength is dominated by ERE. North is to the top and east to the left. Im- 
age is 13" x 73". Note the incipient appearance of the ladder structure that 
is so clear in the HST images, and its sharpening, along with the bright knots 
at the ends of the ladder steps, in Fig. 5 (right image), when the scattered light 
is divided out. 



can actually be combined in the light of the results of a re- 
cently completed investigation (Witt et al. 2006) of the ERE 
in the reflection nebula NGC 7023. This study concluded that 
the ERE results from a two-step process. Step 1 consists of 
the creation of the ERE carrier. This involves the ionization 
(double ionization in case of PAHs) of the ERE carrier pro- 



genitor by photons with energies 10.5 eV < E < 13.6 eV. 
This step can proceed only inside the outflow cavity and its 
walls with a direct line of sight to the central stars, in par- 
ticular to the source of far-ultraviolet photons associated with 
them. Step 2 consists of the pumping of the ERE carrier, fol- 
lowed by fluorescence. The pumping relies on the presence 
of abundant optical/near-UV photons from the luminous AGB 
star HD 44179 located at the center of the RR. 

The evidence for a source of ionizing far-ultraviolet ra- 
diation near the center of the RR, although indirect, is 
strong. The central source in the RR has been identi- 
fied as a binary (Waelkens et al. 1996; Bakkeretal. 1998; 
Men'shchikov et al. 2002) with a orbital perios of 320 days. 
While hidden from direct view at optical wavelengths, the 
center of the RR has been mapped at cm-radio wavelengths 
by Jura et al. (1997) and was found to exhibit a spectrum con- 
sistent with a compact HII region. This is further supported 
by the observations of scattered emission in the Hel recombi- 
nation line at 1.0830 /im by Kelly & Latter (1995), of strong 
hydrogen Balmer alpha line emission (Jura et al. 1997), and a 
host of other permitted and forbidden atomic emission lines 
(Hobbs et al. 2004), all associated with the central region of 
the RR. 

The source of this ionization has remained a subject of some 
controversy. The absence of an observed far-UV excess in the 
spectrum of HD 44179 allowed Men'shchikov et al. (2002) to 
place an upper limit of 100 L on the luminosity of the UV 
source, which, together with information on the mass of the 
secondary star in HD 44179, led them to suggest a hot helium 
white dwarf as the source in question. While it seems clear 
that the photosphere of the Al III primary in HD 44179 can 
be ruled out as the source of the ionizing radiation, Jura et al. 
(1997) discussed the possibility that accretion on either the 
primary or the secondary star in HD 44179 or onto an accre- 
tion disk surrounding the secondary could be the source of 
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Fig. 5 . — Two ratio images, (left panel): 6400 A image divided by 5700 A image. The dark, vertical band in the center is caused by bleeding of the over-exposed 
stellar image, (right panel): 6400 A image divided by 4050 A image. Both images are 73" x 73", and north is to the top and east to the left. The 5700 A and 
4050 A bands are dominated by scattered light, which clearly dominates outside the outflow cones, whereas the ERE dominates in the cavity walls. 



the ionizing radiation. Regardless of whether due to a hot 
white dwarf or accretion, UV continuum radiation in the en- 
ergy range 10.5eV<E< 13.6 eV will escape from the cen- 
tral region preferentially in the direction perpendicular to the 
circumstellar disk, where it can be absorbed by material in the 
walls of the outflow cones to lead to the initiation of the ERE 
carrier. 

While optical/near-UV photons from the Al III primary 
of HD 44179 permeate other parts of the nebula as well, as 
demonstrated by the RR images at wavelengths dominated by 
scattered light, they can excite the ERE only in those places 
where the earlier initiation has created the carriers. Neutral 
PAHs, such as those postulated in the RR disk shadow will 
also absorb these optical/near-UV photons and respond with 
blue fluorescence. If PAHs are responsible for the ERE, the 
initiation constraint found in NGC 7023 suggests that PAHs 
are ionized to the di-cation stage before they luminesce in the 
ERE band. Quantum-chemical calculations of the electronic 
energy level structure of PAH di-cations show the presence of 
strong absorption bands shortward of ~ 6000 A (Witt et al. 
2006). Whether or not PAH di-cations fluoresce in the ERE 
band (> 6000 A) is currently the subject of laboratory inves- 
tigations. 

4.2. Chemical Environments in the Red Rectangle 

Waters et al. (1998) used ISO-SWS spectra from a large 
aperture which included most of the RR nebula (14" x 20") 
for their study of the spectrum of the RR in the mid-IR. These 
spectra showed crystalline silicate emission features and CO2 
absorption features. They also compared broad-band 10 /im 
(N-band) and narrow-band 11.3 /im images, the latter emis- 
sion suggested to arise from a mixture of neutral and ionized 
PAHs (Allamandola et al. 1999; Bakes et al. 2001). These im- 
ages showed that while the 10 /im flux distribution was cen- 
trally peaked, the 11.3 /im image showed an extended distri- 
bution. Bregman et al. (1993) noted that the 11.3 /im image 



showed a narrow-waisted, bipolar shape similar to the ERE 
image. In contrast, their 3.3 /im (arising from predominantly 
neutral PAHs) image was more centrally peaked and lack- 
ing in extended emission in the outflow regions, although the 
3.3 /im is by no means confined to the disk region only. 

Based on the N-band and 11.3 /im images Waters et al. 
(1998) claimed that the disk in the Red Rectangle is made up 
of oxygen-rich material produced in an earlier mass-ejection 
epoch, while the bi-polar outflow lobes are carbon-rich mat- 
ter produced during the more recent outflow. However, the 
spatial resolution of their spectra was insufficient to discrimi- 
nate between the emission from the disk and outflow regions. 
Their claim relies mainly on the argument that the present- 
day outflow is C-rich and that crystalline silicates are unlikely 
to form in an outflow. They concluded, therefore, that the O- 
rich material must be in a stable configuration like the disk. 
Our observations raise questions about such a simple divi- 
sion into two chemically different environments. If the BL 
arises from PAHs (Vijh et al. 2004, 2005) our observations 
demonstrate that such carbon-based molecules are abundantly 
present in the disk as well. Mulas et al. (2006) compute the 
absolute fluxes for the vibrational IR emission from small, 
neutral PAHs carriers proposed by Vijh et al. (2005) and show 
that if anthracene were responsible for the BL, then it would 
almost completely account for the observed flux in the 3.3 /im 
band. Taking into consideration that only 10-20% of the Si 
needs to be crystalline to account for the spectra (Waters et al. 
1998; Kemper et al. 2004) overcomes the argument for all the 
C-rich material to be confined to the outflow cones and the 
disk be completely formed by O-rich material. Our data are 
not inconsistent with the picture of two successive mass ejec- 
tion epochs, a C-rich one following an earlier O-rich one, pro- 
posed by Waters et al. (1998). They do suggest, however, that 
small PAH molecules (if they are the BL carriers) formed dur- 
ing the present C-rich ejection have also entered the space out- 
side the outflow cones. It is in this part of space, especially in 
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the shadow of the circum-binary disk of the RR, where such 
molecules are protected against ionization and are thus able 
to fluoresce in the blue part of the spectrum. 

In summary, the narrow-band morphology of the RR in 
the BL, ERE, near- and mid-IR emission bands gives us fur- 
ther insight into the distribution of chemical environments in 
this system. The disk regions show spectral emission fea- 
tures likely produced by neutral PAHs (3.3 /im emission, BL) 
in addition to crystalline olivine features seen at 11.27 /im 
(Miyata et al. 2004). The outflow regions show emission from 
11.3 /im (Bregman et al. 1993; Miyata et al. 2004) and from 
0.6 /im (the ERE). The former may probably originate from 
neutral and/or singly ionized molecules and the latter from 
doubly ionized PAH species. 

5. CONCLUSIONS 

1. The morphology of the BL and ERE emissions in the 
RR nebula are almost mutually exclusive. The ERE is 
strongly enhanced along the walls and in the interior 
of the bipolar outflow cones of the RR nebula, while 
the BL is strongest in those regions that are located in 
the shadow of the circumstellar disk surrounding the 
central source, HD 44174. 

2. The presence of a circumstellar disk, viewed edge-on, 
divides the space around HD 44179 into two environ- 
ments characterized by radiation fields of greatly dif- 
ferent intensity and hardness. This leads us to suggest 
that the ERE and the BL arise from different ioniza- 
tion stages of the same family of carriers, such as PAH 
molecules. In case of PAH molecules, a large range of 
sizes and species are probably present. The BL pref- 
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erentially detects small, neutral molecules while the 
mid-IR aromatic emission features arise from both neu- 
tral and ionised species. The ERE, if due to PAHs, 
arises from doubly ionized molecules with masses > 
200 amu. 

3. Our observations at 3934 A near the peak of the BL 
band, show an intensity distribution elongated in E-W 
direction, with a maximum at the position of the cir- 
cumstellar disk of HD 44179. If this structure is asso- 
ciated with the outer portions of the disk, its alignment 
differs from that of the inner disk observed in absorp- 
tion at near-IR wavelengths and the disk seen in CO 
emission by 1 1 degrees. 

4. Our observations add to the growing evidence for 
greater chemical complexity of the RR nebula than pre- 
viously proposed. 
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